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Figure  1:  Forces  and  moments  generated  at  the  bound¬ 
aries  of  a  Bemoulli-Euler  beam. 
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1  Introduction 

Derivation  of  the  exact  partial  differential  equations  de¬ 
scribing  motion  of  systems  with  flexible  elements  is  at 
best  tedious  and  at  worst  so  difficult  as  to  drive  the 
analyst  toward  other  pursuits.  A  method  is  presented 
here  which  simplifies  the  process  of  equation  develop¬ 
ment  for  a  wide  variety  of  system  configurations.  The 
key  to  the  method  involves  modeling  each  fiexible  com¬ 
ponent  in  the  system  as  a  single  dynamic  element  which 
captures  the  essential  characteristics  of  the  flexible  mem¬ 
ber.  This  “surrogate  element”,  which  can  be  treated 
substantially  as  if  it  were  a  rigid  body,  stands  in  for  the 
flexible  member  through  most  of  the  derivation.  Any  of 
the  standard  continuum  models  may  be  used  in  the  form 
of  a  surrogate  element  (i.e.  Bemoulli-Euler  beam,  Tim¬ 
oshenko  beam,  etc.)  The  method  also  works  for  systems 
with  rigid  body  rotations.  [Singer  92]  demonstrates  the 
use  of  the  surrogate  element  using  Kane’s  Method. 

First,  the  dynamic  system  is  modeled  with  surro¬ 
gate  elements  substituted  for  flexible  members.  State 
variables  are  then  assigned.  Next,  any  of  the  conven¬ 
tional  methods  for  deriving  system  equations  can  be  ap¬ 
plied.  When  Lagrange’s  method  is  employed,  use  of  sur¬ 
rogate  elements  simplifies  mathematical  manipulations 
and  eliminates  the  need  to  apply  integration  by  parts. 
(Integration  by  parts  is,  generally,  the  most  time  con¬ 
suming  step  in  the  solution  process  for  continuous  sys¬ 
tems.)  The  equation  generation  process  becomes  very 
similar  to  that  for  systems  made  up  of  only  rigid  ele¬ 
ments.  This  representation  and  derivation  strategy  was 
developed  to  facilitate  computer  generation  of  equations 
of  motion  for  flexible  systems.  Its  use  also  provides  ad¬ 
vantages  when  the  equations  are  generated  by  the  ana¬ 
lyst. 

2  The  Surrogate  Element 

A  flexible  component  in  a  dynamical  system  will  typ¬ 
ically  exhibit  the  characteristics  of  both  a  mass  and  a 
spring.  It  will  contain  kinetic  and  potential  energy  and 
will  exert  forces  and  moments  on  adjacent  components 
at  its  boundaries.  If  a  flexible  component  is  beamlike,  it 
is  frequently  represented  as  a  series  of  lumped  elements, 
each  with  specified  degrees  of  freedom.  Bemoulli-Euler 
Beam  elements  and  Timoshenko  Beam  elements  are  two 
common  element  representations.  Each  element  of  the 
series  exerts  a  force  and  a  moment  only  on  its  immediate 
neighbors.  These  forces  and  moments  result  from  local 
deformations  within  the  bodies.  More  specifically,  the 


Figure  2:  This  is  an  example  of  a  surrogate  element. 


load  applied  by  an  element  depends  only  on  the  shape  of 
the  element  at  the  point  of  application.  (Of  course,  from 
a  slightly  different  perspective,  the  shape  is  dependent 
only  on  the  load.)  An  arbitrary  number  of  similar  ele¬ 
ments  can  be  chosen  to  represent  a  flexible  component. 

Just  as  the  load  apphed  by  an  element  onto  its  neigh¬ 
boring  element  depends  only  on  its  shape  at  the  bound¬ 
ary,  the  forces  and  moments  apphed  by  the  entire  flexible 
component  onto  adjacent  lumped  or  flexible  components 
are  functions  only  of  the  flexible  component’s  shape  (spa¬ 
tial  derivatives)  at  its  boundaries.  For  a  Bemoulh-Euler 
Beam  of  length  I,  pictured  in  Figure  1,  these  forces  and 
moments  are  [Crandall  82] 
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The  sign  convention  chosen  for  these  equations  is  as  set 
by  the  directions  of  the  vectors  in  Figure  1. 

The  surrogate  element  is  a  model  or  representation  of 
a  flexible  component  which  can  be  lased  instead  of  a  tra¬ 
ditional  continuous  representation  during  equation  gen¬ 
eration.  The  surrogate  element  exhibits  the  essential  dy¬ 
namic  characteristics  of  the  component.  It  consists  of  a 
single  elemental  particle  of  the  continuum  (called  a  “rep¬ 
resentative  particle”)  and  a  spine  defining  the  element’s 
shape.  The  element  can  exert  forces  and  moments  at 
the  boundaries  that  connect  it  to  the  rest  of  the  system. 
The  surrogate  element  contributes  a  kinetic  energy  term 
because  of  motion  of  its  elemental  mass,  and  a  potential 
energy  term  due  to  energy  storage  at  the  location  of  the 
representative  particle.  Because  every  internal  particle 
of  the  continuum  is  topologically  identical,  the  response 
of  one  particle  can  be  used  to  represent  that  of  the  con¬ 
tinuum.  The  forces  and  moments  at  the  boundaries  en¬ 
force  the  connection  of  the  surrogate  element  to  the  rest 
of  the  system.  (Appendix  A  justifies  this  description  of 
the  surrogate  element.) 

Generalized  coordinates  must  be  assigned  to  the  sur¬ 
rogate  element  at  each  of  the  two  boundary  points  and 


at  the  representative  particle  (Table  1  shows  coordinate 
selection  for  a  variety  of  surrogate  elements.)  Frequently 
the  generalized  coordinates  which  locate  the  ends  of  the 
surrogate  element  will  also  locate  the  bodies  to  which 
the  corresponding  beam  is  rigidly  attached. 

To  satisfy  constraints  set  by  Lagrange’s  method,  the 
generalized  coordinates  locating  the  representative  par¬ 
ticle  must  be  referenced  to  a  coordinate  system  which 
is  geometrically  compatible  with  that  of  the  bodies  to 
which  the  surrogate  element  is  connected.  For  example, 
in  Figure  2: 


turn,  commimicate  with  the  surrogate  elements  through 
the  enforcement  of  boimdary  position  and  orientation 
compatibility.  The  loads  applied  by  adjacent  compo¬ 
nents  on  the  surrogate  element  are  captured  implicitly 
by  these  geometric  constraints.  Consequently,  they  are 
not  included  explicitly  in  the  derivation  (see  Appendix 

A). 

3  An  Extended  Form  of  Lagrange’s 
Equations 
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The  condition  of  geometric  compatibiUty  is  discussed  in 
Appendix  B.  Note  that  the  surrogate  element  in  Figure  2 
represents  a  Bemoulli-Euler  beam  and  so  the  position  of 
the  representative  point  can  be  set  with  a  single  gen¬ 
eralized  coordinate  [Crandall  82].  If  this  were  a  surro¬ 
gate  element  for  a  Timoshenko  beam,  two  generalized 
coordinates  would  be  required  to  set  both  position  and 
orientation  at  that  point. 

While  the  Bernoulli-Euler  beam  has  only  one  coor¬ 
dinate  for  describing  the  motion  of  each  beam  element, 
two  coordinates  are  required  at  each  boundary  to  specify 
position  and  orientation  of  the  end  of  the  beam.  This 
fact  is  consistent  with  the  continuous  derivation  of  the 
Bernoulli-Euler  beam  model  [Crandall  82]  and  is  not  a 
product  of  the  surrogate  element  definition.  (The  surro¬ 
gate  element  concept  simply  highlights  this  anomaly.) 

For  the  case  of  the  Bemoulli-Euler  Beam,  the  surro¬ 
gate  element  has  kinetic  coenergy,  potential  energy,  and 
Lagrangian  given  by 


r  -  (6) 


L,  =  T*-V  (8) 


with  terms  as  defined  at  the  bottom  of  Table  1. 

Note  that  the  kinetic  energy  term  contains  the  ab¬ 
solute  velocity  of  the  surrogate  element’s  representative 
particle,  Vabs-  If  the  neutral  axis  of  the  beam  is  sta¬ 
tionery,  then  Vabs  —  ^  •  However,  if  the  beam  is  rotating 
or  translating,  the  absolute  velocity  of  the  representa¬ 
tive  particle  must  be  determined  (see  Example  1).  This 
is  consistent  with  the  specification  of  the  kinetic  energy 
of  a  lumped  mass.  Table  1  includes  expressions  for  ki¬ 
netic  and  potential  energy  for  several  types  of  surrogate 
elements. 

In  the  derivation  strategy  developed  below,  surrogate 
elements  are  assumed  to  communicate  with  adjacent 
lumped  components  by  applying  loads  to  them.  Table  1 
includes  expressions  for  these  loads  as  they  should  ap¬ 
pear  in  the  equations.  The  lumped  components  will,  in 


2 


The  common  expression  of  Lagrange’s  equations: 


dt  ~ 


i=l,2,...,n  (9) 


is  valid  for  Lagrangians  of  the  form, 
i(6,  Here  Ej  (for  j  =  1, 2, . . . , n) 

are  the  generafized  forces  [Crandall  82].  When  a  system 
model  includes  surrogate  elements,  the  potential  energy 
terms  become  functions  of  the  spatial  derivatives  of  the 
^i',  ^  and  Appendix  A  demonstrates  that  a  special 
form  of  Lagrange’s  equations, 


(  dL,  \ 


>Ej  j  =  l,2,...,n  (10) 


can  be  derived  for  systems  whose  models  include  rigid 
and  surrogate  elements.  For  this  equation,  X,  is  the  La¬ 
grangian  formulated  to  include  surrogate  element  terms 
and  the  ’Ej  (for  j  =  1,2,  ...,n)  are  the  generalized 
forces  which  include  the  loads  applied  by  the  surro¬ 
gate  elements  onto  adjacent  components.  The  first  two 
terms  of  Equation  10  are  those  of  Equation  9.  The 
last  two  terms  are  those  required  when  spatial  deriva¬ 
tives  appear  in  the  Lagrangian.  This  result  is  shown 
by  [Lanczos  70]  and  [Meirovitch  90]  and  is  also  derived 
for  use  in  Lagrangian  density  functions  by  Low  and 
Vidyasagar  [Low  88]. 


4  Derivation  Strategy 

The  steps  in  a  strategy  for  using  sm^ogate  elements  to 
derive  exact  partial  differential  equations  of  motion  for 
systems  with  flexible  components  will  be  described  here. 
The  description  assumes  the  use  of  Lagrange’s  equation 
and  is  followed  by  a  sample  derivation  chosen  to  illustrate 
application  of  the  method. 

1.  Model  the  system  as  a  collection  of  rigid  and  flexi¬ 
ble  elements. 

2.  Temporarily  substitute  for  each  flexible  component 
a  surrogate  element. 

3.  Assign  the  required  generalized  coordinates  to  the 
surrogate  elements.  If  the  system  configuration 
precludes  motion  in  any  direction,  a  generalized 
coordinate  in  that  direction  need  not  be  included. 


4.  Assign  any  additional  generalized  coordinates  re¬ 
quired  to  specify  the  states  of  the  system’s  rigid 
elements.  ^ 

5.  List  the  kinetic  and  potential  energies  of  all  surro¬ 
gate  elements  as  presented  in  Table  1. 

6.  List  the  kinetic  and  potential  energies  of  the 
lumped  elements. 

7.  From  these  terms  for  kinetic  and  potential  energy, 
form  a  Lagrangian. 

8.  List  the  generalized  forces.  The  portions  of  the  gen¬ 
eralized  forces  applied  by  the  surrogate  elements 
to  adjacent  components  are  given  in  Table  1. 

9.  Apply  the  extended  form  of  Lagrange’s  equations 
to  the  Lagrangian  for  each  state  variable. 

10.  For  those  resulting  equations  for  which  the  variable 
specifying  the  width  of  the  representative  particle 
is  a  conunon  factor,  divide  through  by  it.  In  the 
example  below,  this  variable  is  Ax.  Note  that  these 
equations  will  be  those  for  generalized  coordinates 
representing  displacements  of  representative  parti¬ 
cles  from  their  neutral  positions. 

11.  Integrate  along  the  spine  of  each  surrogate  element 
those  terms  of  the  equations  which  still  contain  Ax, 
the  variable  specifying  the  width  of  the  representa¬ 
tive  particle  of  the  surrogate  element.  Performing 
this  step  removes  the  surrogate  elements  from  the 
analysis  and  3aelds  the  exact,  partial  differential 
equations  of  motion  and  the  boundary  conditions 
for  the  system  as  modeled  in  step  1. 

5  Example  One 

A  flexible  beam  is  fixed  to  a  moving  mass,  mj,  as  shown 
in  Figure  3a.  A  mass,  mo,  is  coimected  to  ground  with  a 
spring,  mb  is  constrained  to  move  only  vertically  and  is 
connected  with  a  spring  and  a  damper  to  mo-  The  mass, 
me  is  connected  to  groimd  through  a  third  spring.  A 
distributed  force,  /(x),  is  applied  to  the  beam.  System 
parameters  are  as  presented  in  the  figure. 

Figiure  3b  shows  the  system  with  a  surrogate  element 
substituted  for  the  flexible  beam.  Because  the  left  end  of 
the  surrogate  element  spine  cannot  rotate,  a  generalized 
coordinate  is  assigned  for  translation  {q^)  but  not  rota¬ 
tion.  Coordinates  and  qs  represent  the  displacement 
and  rotation  of  the  right  end  of  the  beam  respectively. 
Coordinate  qi  locates  mass,  m^.  Note  that  gs,  and  qi 
are  all  referenced  to  the  same  coordinate  system;  they 
all  have  been  defined  relative  to  mass,  ma- 

The  independent  coordinates  are  chosen  as; 

^i  =  qi  i=l,2,...,5  (11) 


a)  Model  of  system  with  both  rigid  body  and  continuous 
elements. 


b)  Model  of  system  with  surrogate  element  and  assigned 
coordinates. 


^Be  certain  that  the  coordinates  for  the  parts  of  the  rigid 
system  that  connect  to  each  surrogate  element  are  defined  to 
be  compatible  with  the  coordinates  for  the  corresponding  con¬ 
tinuum.  Appendix  B  elaborates  on  this  requirement. 


Figure  3:  System  configuration  for  Example  1. 
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The  kinetic  coenergy,  potential  energy,  and  La- 
grangian  are  expressed  as: 


2-^191^  +  2-^292^  +  2^3  (91  +  94)^  (13) 


L,  =  T*-V  (14) 

Note  that  V  is  now  a  function  of  ^  (or  ^). 

The  new  generalized  forces,  ‘Ei,  as  defined  in  Sec¬ 
tion  3  are  determined  for  each  degree-of-freedom: 


'El  =  /(x)Arr  (15) 

'E2  =  (16) 

'Es  =  fix)Ax  (17) 

'H4  =  Fi  (18) 

'Es  =  M2  (19) 


For  this  example  the  terms  to  be  substituted  into  La¬ 
grange’s  equations  (10)  for  ^1  =  91  are 
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Applying  Equation  10  yields 
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Dividing  Equation  27  by  Ax  (Step  10)  yields 
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6  Conclusions 
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A  method  has  been  proposed  for  deriving  the  partial 
differential  equations  of  motion  for  systems  of  lumped 
and  flexible  components.  When  used  in  conjunction  with 
Lagrange’s  method,  the  process  of  equation  generation 
is  significantly  simplified.  Particularly,  integration  by 
parts  is  no  longer  required.  Both  system  equations  and 
expressions  of  boundary  conditions  result  directly  fi:om 
application  of  this  method.  The  method  contains  the 
same  simple  steps  used  to  derive  the  equations  of  motion 
for  lumped  parameter  systems. 


7  Appendix  A 


Terms  in  Equation  25  can  then  be  integrated  over  the 
length  of  the  beam  (Step  11).  For  example. 
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where  rribeam  is  the  total  mass  of  the  beam. 

The  final  set  of  partial  differential  equations  (after 
Step  11)  is: 
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=  -b^ 


f{x)dx 

dqz 

^qz 


(32) 


dt 

dx 


El 


9x2 


+ 


^  9Z2  +  dt^  ) 

^  \-.r^qz 


9x2 


9x2 


=  f{x) 


(33) 

x=0 


(34) 


(d^q\  ,  A4' 

*=  V  dt^  d^  . 


In  this  section  it  will  be  shown  that  the  application  of 
the  surrogate  element  method  yields  correct  equations 
of  motion  and  boimdary  conditions  when  apphed  to  a 
dynamic  system  of  rigid  and  flexible  elements.  First, 
equations  will  be  derived  for  a  general  system  using  the 
variational  formulation.  Then  the  surrogate  element  for¬ 
mulation  will  be  employed,  yielding  identical  equations. 

A  Lagrangian  for  a  system  of  rigid  and  flexible  com¬ 
ponents  can  be  express^  as, 

n 

L  =  ^  ^  Lflexiblej  d"  Lrigid  (3'^) 

J=1 

where  LflexiWej  contains  all  of  the  kinetic  and  potential 
energy  terms  for  the  jith  of  n  continuous  elements,  and 
Lrf^d  contains  the  rest  of  the  terms,  those  not  dependent 
on  spatial  orientation  within  an  element.  LBadhiej  is  the 
sum  of  the  energy  terms  for  the  elemental  particles  in 
the  jth  continuum,  so  the  Lagrangian  can  be  written 

n 

fjAzj  -J-  (38) 

j=l  Zj 

where  the  Lij  is  a  Lagrangian  density  function 
(see  [Low  88]  and  [Meirovitch  90])  and  Zj  is  a  variable 
expressing  the  position  of  each  element  along  the  jth 
continuum. 

The  components  of  the  Lagrangian  for  the  system  are 
expressed  in  functional  form  as 

Lf j  =  Lf j  (rj ,  Tj ,  &j {zj),  Cj {zj), 

ei(2j),eJ(2j),t)  (39) 

•Lrigid  —  -^ripd  (*"«>  Xi,  ej(2:jg),  ej(^^jj), 

where  the  n  are  state  variables  defining  the  positions  of 
the  m  rigid  elements  (for  i  =  1, 2, . . . ,  m);  the  ej{zj)  are 
coordinates  which  define  deflections  within  the  flexible 
elements;  Zj^  and  Zj^  axe  the  locations  of  the  endpoints 
of  the  jth  continuous  element;  and  ej{zj^)  and  ej{zj^) 


Table  1:  This  is  a  table  of  example  surrogate  elements. 

6 


are  the  coordinates  at  those  endpoints.  Differentiation 
in  time  is  denoted  with  a  dot  {e.j{zj)  =  dej{zj)/dt)-,  spa¬ 
tial  differentiation  is  denoted  with  a  prime  (e'(zj)  = 
dej{zj)ldzj).  The  time  derivative  of  e'j{zj)  is  denoted 
as  ^{zj). 

The  L[j  functions  contain  the  terms  in  L  which  are  to 
be  spatially  integrated,  while  the  Lripd  function  includes 
the  remaining  terms.  The  Lrigid  function  contains  con¬ 
tinuous  coordinates  specified  at  the  boundaries  of  the 
continuimi  (the  ej{zj^)  and  ej{zj^)).  For  a  system  of 
only  rigid  elements,  a  single  variable  would  be  selected 
to  specify  the  motion  of  each  rigid  component.  However, 
for  mixed  systems,  coordinates  assigned  to  continuous  el¬ 
ements  specify  the  motion  of  some  rigid  components  of 
the  system  so  additional  variables  are  not  needed.  For 
example,  in  Figure  3,  ©(0)  =  ga-  For  the  variational 
approach,  53(0)  appears  in  Ltigid-  With  the  smrogate 
element  approach,  ©  would  be  used  instead. 


7.1  The  Variational  Approach 

The  summation  of  the  energy  terms  in  Equation  38  be¬ 
comes  an  integral  when  the  continuous  case  is  consid¬ 
ered, 

L  = 

The  Variational  Indicator  may  be  formed  from  the 
Lagrangian,  L  and  the  work,  W,  done  by  generalized 
forces; 


u 


IjfjdZj  -|- 


(41) 


ot2 

j  ' 


SI  {L  +  W)dt  = 

y  SLtj{ri,ri,ej{zj),ij{zj), 

A  . 

"h  /  ^-l^ripd  (^ii 
Jti 


pt2 

+  /  <  y^Srigidj^rj-ny^  f  Efjdzj  6ej{zj) 

J  tx  i  i  ’ 


)dt  (42) 


The  Ejigidj  are  the  generalized  forces  in  the  Sri  directions 
(associated  with  the  rigid  body  coordinates).  The  Ef^ 
are  the  generalized  forces  in  the  Sej{zj)  directions  (asso¬ 
ciated  with  the  flexible  coordinates).  The  Ebdj,  EbcZi) 
Ebc3j,  and  Ebc4j  are  the  generalized  forces  in  the  Sej{zj^, 
Sej{zjJ,  Se'j{zj^),  and  Se'j{zjJ  directions  respectively. 


These  generalized  forces  are  derived  using  the  Principle 
of  Virtual  Work,  and  therefore  do  not  contain  internal 
force  terms.  The  subscript  Ebci  •  ■  •  Ebc4  refer  to  the  E’s 
associated  with  the  four  boimdary  coordinates  of  each 
flexible  element. 

The  variation  operator  is  apphed  to  each  term. 


E' 


+ 


+ 


'  '  dh 

j  Jzj 

t-  /  /  Etjdzj  Sej{zj) 

J  Zj  ^  ^  J  Zj 


I  ^-^rigid  c  /  ^  ^J'rigid 

+“bc3j^Cjf(^j{,)  +  “bc4j^^j(^Jo) 


Next,  the  variation  is  integrated  by  parts.  The  terms 
in  the  variation  of  time  derivatives  are  integrated  first  in 
time.  The  terms  in  the  variation  of  spatial  derivatives 
are  integrated  first  in  Zj, 


J  t\  V  t\  ^  j 


+ 


Z—^  /  dn 

j  ^ 


•i  J  Zi 
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Sri 


The  next  set  of  equations  is  derived  from  setting  the 
coefficient  of  each  6ej{zj)  to  zero.  The  summation  over  j 
disappears  because  each  term  of  the  summation  contains 
a  different  6ej{zj)  for  any  given  j,  and  therefore,  becomes 
a  separate  equation.  For  the  integral  in  Zj  to  be  zero  for 
arbitrary  variations,  the  integrand  must  be  zero  yielding 


d  (  dLrj  \ 
dt  [dej(zj)j 


de'.{zj) 


dej{zj) 


dzj^  \de'!{zj) )  ’ 

j=l,2,...,n  .  (46) 

The  last  four  sets  of  equations  are  boundary  condi¬ 
tions  at  the  edges  of  the  continuum  corresponding  to 
Note  that  the  summation  over  j  is 
dropped  for  the  same  reason  as  in  the  previous  equation. 


dt  \dti{z^,))  ~ 


=  Hbcij  j  =  1,2,  ...,n  (47) 


f  ^-^rigid  \ 


dLfj 


—  ^bcZj 


.,n  (48) 


Note  that  the  summation  in  i  is  for  each  of  the  rigid 
body  equation  terms;  the  summation  over  j  is  for  the 
terms  due  to  the  flexible  components.  Now  the  equations 
of  motion  can  be  extracted  from  Equation  44.  The  first 
set  of  equations  arises  from  setting  the  coefficient  of  each 
6ri  to  zero. 


—  (  ^^rigid  \ 

dt  \^ae;.(0,^)  j 

dL(j 

de'd^ib) 


^■^rigid 

HbcSj  J  =  1,2, . . .  ,n  (49) 


(  ^4igid  \ 
*  ) 

I  dU,  _ 


^^rigid 

d^d^ja) 

•^bc4j  3  ~  1,  2,  .  .  .  ,  7Z  (50) 


In  total  there  are  m  rigid  body  equations,  n  flexible  com¬ 
ponent  equations,  and  4n  boundary  equations. 


7.2  The  Surrogate  Element  Method 


The  equation  derivation  presented  above  is  now  per¬ 
formed  with  the  same  notation  but  using  the  surrogate 
element  method  outlined  in  this  paper.  The  Lagrangian 
for  this  case  is 

n 

^  ^  ]  djfj^Zj  -f-  Lrigid  (51) 

j=l  Zj 

with  terms  as  defined  in  Equations  39  and  40.  The  La¬ 
grangian  is  formed  from  the  equivalent  rigid-particle  sur¬ 
rogate  element  model.  Equation  10  is  applied  to  the 
Lagrangian, 


d  (  dL,\ 
^\diJ  d^,^dzj  (^a(||)J 


(52) 


-1 

i 


The  first  m  of  the  correspond  to  coordinates  n  above; 
the  next  n  of  the  correspond  to  coordinates  ej(zj); 
the  last  4n  of  the  correspond  to  coordinates 

respectively.  Therefore,  p  = 
m  +  5n,  and  j  =  k  —  m  so  that  the  zj  refer  to  the  spatial 
variables  along  the  continuum  corresponding  to  ^k-  The 
®Efc  are  the  generalized  forces  defined  in  terms  of  the 
generalized  forces  from  the  variational  approach  given 

by 


“rigidfe 

m) 


"bcl  (/s— m— n) 
+'^tablefc 


“ bc2  (A:— m— 2n) 

"I"  *^^tabl€A: 


bc3  (A;— m— 3n) 
+“tabl«Jb 


“bc4(A;— m— 4n) 
+5tablei 


for  k  = 

1,2,. ..,m 
for  k  = 
m+  1, 

m  +  2, . . .  ,m  +  n 
for  k  = 
m  +  n+1, 
m  +  n  +  2, 

. . . ,  m  +  2n 
for  k  = 
m  +  2n  +  1, 

T7l  272  +  2, 

. . . ,  m  +  3n 
for  k  = 
m  +  372+  1, 

772  +  372  +  2, 

.  .  .  ,  7T2  +  472 
for  k  = 

772  +  472+  1, 

772  +  472  +  2, 

.  .  . ,  772  +  572 


(53) 


where  the  Stable*  are  obtained  from  Table  1  and  can  only 
be  non-zero  for  coordinates  at  the  boundaries  of  the  con- 
tinua. 

The  correspond  to  three  different  variable  types; 
the  rigid  body  coordinates,  the  continuous  coordinates, 
and  the  continuous  coordinates  at  the  boundaries.  The 
evaluation  of  Equation  52  is  treated  separately  for  each 
case.  For  A;  =  1, 2, . . . ,  m.  Equation  52  becomes, 

\  \  ''  A  -  ^^rigid 


j  Azj 

d  f  \  _ 

dt  \  dvi  J 


The  third  and  fourth  terms  in  Equation  52  are  zero  be¬ 
cause  the  ^k  are  rigid  coordinates. 

For  k  =  772+1,  772+2, . . . ,  772+72,  Equation  52  becomes. 


dLfj  d  f  dLfj  N 

dej{zj)  dtydcjizj)} 

'  dLfj 


d 


dzj  \de'j{zj) 


A 


dLfj 

de'!{zj) 


Azj 


=  Ef  jAzj  j  =  l,2,...,n 


(55) 


becomes, 

^  __  ^-^rigid  _ 

dt  \dej{zj^))  dcjizj^) 

Hbcij  +  Pij 

d  /  dLjjpi  \  dL^^i  _ 

\dej[zjj)  dejizjj 

“bc2j  "I"  L^j 


dL 


'rigid 


9e'j{zjt) 


dLg^i  _ 

de'jizjb) 

HbcSj  +  Mij 
^Trigjd 


d^jizjj 
Hbc4j  +  M2j 


j  =  l,2,...,n  (56) 


i  =  l,2,...,72  (57) 


j  =  l,2,...,72  (58) 


j  =  1,  2,  .  . .  ,  72  (5^) 


where  the  Fij,  F2j,  Mij,  and  M2j  are  obtained  from 
Table  1.  The  Table  defines  these  values  as  the  forces  and 
moments  acting  from  the  continuum  onto  the  adjacent 
components.  The  values  used  in  creating  the  table  are 
given  by 


Fi  = 
F2  = 
Ml  = 
M2  = 


dL 


fi 


dL 


fi 


de'jizjf,)  dzj  yde'({zj^)^ 


dLfj 

9e'j{zjJ 

dLfj 

de'!{zjk) 

dLfj 

de'jizjj 


_d_  (  dLfj 
dzj  [de'!(zjj 


(60) 

(61) 

(62) 

(63) 


Following  the  procedure  presented  in  Section  4,  Azj  is 
canceled  from  every  term  of  Equation  55  to  yield  Equa¬ 
tion  46.  Next,  all  terms  in  Azj  that  remain  in  any  of  the 
equations  are  integrated  along  the  Zj.  These  terms  are 
all  foimd  in  Equation  54.  The  integration  yields  Equa¬ 
tion  45.  Finally,  Equations  56-59  can  be  rewritten  with 
substitutions  from  Table  1  for  the  Fi,  F2,  Mi,  and  M2 
to  obtain  Equations  47-50. 

8  Appendix  B 

The  surrogate  element  derivation  above  assumed  that 
the  coordinate  describing  the  location  of  the  representa¬ 
tive  particle  in  a  continuum  be  geometrically  compatible 
with  the  coordinates  for  adjacent  elements  at  the  bound¬ 
aries  of  that  continuum.  Geometric  compatibility  means 
that  the  variation  of  the  surrogate  element  coordinates 
and  the  variation  of  the  boundary  coordinates  be  equal. 
Consider  the  example  shown  in  Figiue  4.  Coordinate  r2 
at  X  =  0  equals  ri,  and  r2  at  x  =  Z  equals  rg.  However, 
92  at  X  =  0  is  equal  to  zero,  and  93  equals  the  sum  of  91 
plus  92  at  X  =  Z.  The  method  outlined  above  assumes 
that  the  consistent  coordinate  set  (n,  r2,  rs)  be  chosen. 
The  requirement  on  the  coordinate  selection, 


For  k  =  772+72+1,772+72+2, . . .  ,772+572,  Equation  52 


(64) 


I 


t 


Figure  4;  Example  using  coordinates  with  different  refer¬ 
ence  frames.  The  ri,  r2,  frame  is  referenced  to  ground. 
The  qi  and  qs  are  referenced  to  ground  while  q2  is  refer¬ 
enced  to  the  left  hand  mass. 


and 

=  Sn  (65) 

is  satisfied.  The  force  terms  on  the  end  of  the  contin¬ 
uum  are  multiplied  by  the  variations  of  the  coordinates 
at  those  boundaries  (^r2|^_o  and  Sr2\^^i  in  this  exam¬ 
ple).  The  force  terms  acting  on  the  rigid  components 
are  multiplied  by  the  variations  of  the  rigid  component 
coordinates  (Sri  and  Sr^  in  this  example).  The  method 
outlined  above  implicitly  equates  these  two  sets  of  forces. 
However,  if  the  two  sets  of  variations  are  not  equal  for 
the  chosen  coordinate  system,  then  the  equating  of  these 
force  terms  can  not  be  justified. 

The  requirement  discussed  above  can  be  relaxed  by 
realizing  that  the  force  terms  can  be  corrected  when  us¬ 
ing  a  selection  of  coordinates  that  is  not  geometrically 
compatible.  In  the  example  above  using  the  qj,  q2,  93 
system,  ^q2li=o  =  0.  The  force  term  Fi,  which  would 
be  applied  to  the  mass  on  the  left  when  using  a  geo¬ 
metrically  compatible  coordinate  system,  is  not  added 
to  ®5i.  Also,  =  Sq3  -  Sqi.  Therefore,  for  La¬ 

grange’s  method  to  provide  correct  equations  of  motion, 
F2,  the  force  at  the  right  end  of  the  beam  must  be  added 
to  the  generalized  force  term  for  the  93  equation  (^Es), 
and  subtracted  from  the  generalized  force  term  for  the 
9i  equation  (®Ei).  This  modification  is  convenient  when 
applying  the  surrogate  element  method  to  more  complex 
systems  such  as  those  with  rotating  flexible  components. 
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